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Abstract: 

First numerical results of the Monte Carlo generator RacoonWW for e + e _ — > 
WW — > 4f (+7) in the electroweak Standard Model are presented. This event generator 
is the first one that includes 0{a) electroweak radiative corrections in the double-pole 
approximation completely. We briefly describe the strategy of the calculation and give 
numerical results for total cross sections, including CC03, and various distributions. 
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At present, the focus of Standard-Model (SM) tests lies on W-boson-pair production 
at LEP2, i.e. on the process e + e~ — > WW — > 4f. In order to match the experimental 
accuracy of roughly 1%, theoretical predictions for the cross sections of e + e~ — > 4/ with a 
precision at or below the per-cent level are needed. This requires to include the complete 
set of lowest-order diagrams for e + e~ — > 4/ and the 0(a) corrections to the W-pair 
production channels e + e~ — > WW — > 4/. Since the impact of electroweak corrections 
grows with increasing energy, this task is even more important for future linear colliders 
with higher energy and luminosity. 

In all regions of phase space where W-pair production dominates the cross section 
for e + e~ — > 4f, an expansion of the matrix element about the poles of the resonant 
W propagators provides a reasonable approximation for the radiative corrections. Ne- 
glecting corrections to the non-doubly-resonant contributions leads to uncertainties of the 
order a/it x r w /M w x l°g(. . .) ~ 0.1% with respect to the leading lowest-order contri- 
butions, where the logarithm indicates possible logarithmic enhancements. Thus, this 
so-called double-pole approximation (DPA) should be sufficient for an accuracy of 0.5% 
for observables that are dominated by doubly-resonant diagrams. Moreover, the DPA 
provides a gauge-invariant answer and allows us to use the existing results for on-shell 
W-pair production [ [TJ [| and W-boson decay [ ^ |J, as far as the virtual corrections are 
concerned. 

The DPA for 0(a) corrections to pair production of unstable particles has already 
been used in the literature. A possible strategy has been proposed in Ref. [ [J. A Monte 
Carlo generator including the corrections to the W-pair production subprocess and the 
leading- logarithmic corrections to the W-boson decays has been constructed Refs. [ |6], |7j , 
but non-factorizable corrections and W-spin correlations have been neglected there. A first 
complete calculation of the 0(a) corrections for off-shell W-pair production, including a 
numerical study of leptonic final states, was presented in Ref. [ || using a semi-analytical 
approach, which is, however, only applicable to ideal theoretical situations. 

In this paper we present the first complete calculation of the O(o) corrections for off- 
shell W-pair production in DPA that has been implemented in a Monte Carlo generator, 
which is called RacoonWW. This generator includes the complete lowest-order matrix 
elements for e + e~ — > 4/ for any four-fermion final state. For the virtual corrections a 
DPA is used without any additional approximations. In particular, the exact four-fermion 
phase space is used throughout. The virtual corrections consist of factorizable and non- 
factorizable contributions. The former are the ones that are associated to either W-pair 
production or W-boson decay; the results of Refs. [ [l], [| are used in this part. The latter 
comprise all corrections in which the subprocesses production and decays do not proceed 
independently. Up to some simple supplements, the virtual non-factorizable corrections 
can be read off from the literature [ || |10| ; we made use of the results of Ref. [ |l(J . The real 
bremsstrahlung corrections are based on the full matrix-element calculation for e + e~ — > 



4/7 described in Ref. [ [D] . More precisely, the minimal gauge- invariant subset including 
all doubly-resonant contributions of the processes e + e~ — > WW — > 4/7, i.e. the photon 
radiation from the CC11 subset, are included. By using the exact matrix elements for 
the real radiation, we avoid problems in defining a DPA for semi-soft photons (E 1 ~ r\y) 
and, moreover, can include the leading logarithmic corrections to non-doubly-resonant 
diagrams (background diagrams) exactly. The real corrections and the virtual corrections 
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are matched in such a way that all infrared singularities cancel exactly. The initial-state 
collinear singularities are regularized by retaining a finite electron mass and factorized into 
lowest-order matrix element and splitting functions. The collinear singularities connected 
to final-state radiation are treated inclusively, i.e. photons within collinear cones around 
the final-state fermions are integrated over, so that no logarithmic final-state fermion mass 
dependence survives. All contributions have been implemented in two programs, one of 
which uses the subtraction method described in Ref. [ |T3], the other one uses phase-space 
slicing. All parts of the calculations have been performed in two independent ways. A 
detailed description of the calculation and the Monte Carlo generator RacoonWW will 



be published elsewhere [ 13 



For the numerical results we used the fixed- width scheme and the following parameters: 



G> — 


1.16637 x KT 5 GeV~ 2 , 


a = 


1/137.0359895, 


M w = 


80.35 GeV, 


r w = 


2.08699 . . . GeV, 


M z = 


91.1867 GeV, 


r z = 


2.49471 GeV, 


m t = 


174.17 GeV, 


M H = 


150 GeV, 


m e = 


510.99907 keV. 







The weak mixing angle is fixed by c w = Mw/Mz, = 1 — c^. These parameters are 
over-complete but self-consistent. Instead of a we use to parametrize the lowest-order 
matrix element, i.e. we use the effective coupling 

«G M = (2J 

71 

in the lowest-order matrix element. This parameterization has the advantage that all 
higher-order contributions associated with the running of the electromagnetic coupling 
and the leading universal two-loop m t -dependent corrections are correctly taken into ac- 
count. In the relative 0(a) corrections, on the other hand, we use a, since in the real 
corrections the scale of the real photon is zero. The W-boson width given above is cal- 
culated including the electroweak and QCD one-loop corrections with a s = 0.119. We 
do not include QCD corrections to the process e + e~ — > WW — > 4/, and initial-state 
radiation is only taken into account in 0(a). 

In Tables [[J |2], and [3] we present numbers for total cross sections without any cuts, 
for centre-of-mass (CM) energies 184, 189, and 200 GeV, respectively, based on 20 million 
events. In particular, we give the CC03 cross sections, i.e. the cross sections resulting 
from the signal diagrams only (defined in the 't Hooft-Feynman gauge). We also give 
numbers resulting from the complete set of diagrams for those final states where this is 
possible without cuts, i.e. the CC11 class of processes. In the considered cases, the effects 
of the background diagrams are below 0.2%. Note that we treat the external fermions as 
massless. Therefore, the cross sections for processes not in CC11 class become singular 
if no cuts are imposed for final-state electrons collinear to the beams and for virtual 
photons splitting into // pairs with small invariant masses. The shown corrections are 
typically -13%, -12%, and -9% at 184, 189, and 200 GeV, respectively. The numbers 
in parentheses are estimates of the Monte Carlo integration errors. The numbers for the 
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final state 


CC03 Born 


full Born 


CC03 corrected 


full corrected 




210.26(5) 


210.57(5) 


182.21(11) 


182.53(11) 


z/^/i+du 


630.8(2) 


631.8(2) 


546.5(3) 


547.4(3) 


udsc 


1892.3(5) 


1895.4(5) 


1637.7(7) 


1640.7(7) 


total 


17031(3) 




14749(5) 





Table 1: Total cross sections in fb for e + e — > WW — > 4/ without cuts for various final 
states at 184 GeV 



final state CC03 Born full Born CC03 corrected full corrected 





216.07(6) 


216.39(6) 


190.64(11) 


190.96(11) 


z/^/i+du 


648.2(2) 


649.2(2) 


571.6(3) 


572.6(3) 


udsc 


1944.7(5) 


1947.7(5) 


1713.4(8) 


1716.4(8) 


total 


17501(3) 




15429(5) 





Table 2: Total cross sections in fb for e + e — > WW — > 4/ without cuts for various final 
states at 189 GeV 



final state 


CC03 Born 


full Born 


CC03 corrected 


full corrected 




219.82(6) 


220.06(6) 


199.74(12) 


199.98(12) 


z/^/i + du 


659.5(2) 


660.2(2) 


598.9(3) 


599.7(3) 


udsc 


1978.4(5) 


1980.8(5) 


1795.0(9) 


1797.4(9) 


total 


17806(3) 




16164(5) 





Table 3: Total cross sections in fb for e + e — > WW — > 4/ without cuts for various final 
states at 200 GeV 



total W-pair production cross section in the last row of the tables directly result from the 
other rows by multiplying these with the number of equivalent channels and adding them 
up. 

Next, we study various angular distributions. Here we restrict ourselves to the z/^/z+du 
final state and to yfs = 200 GeV. Because of our treatment of mass singularities (see 



Ref. [ 13] for details) it is necessary to combine photons that are collinear to the incoming 
or outgoing fermions appropriately with these fermions in order to obtain well-defined 
finite distributions. To this end we introduce the following recombination and cut proce- 
dure which proceeds in three steps: 

1. All photons within a cone of 5 degrees around the beams are treated as invisible, 
i.e. their momenta are disregarded when calculating angles, energies, and invariant 
masses. 
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2. Next, the invariant masses of the photon with each of the charged final-state fermions 
are calculated. If the smallest one is smaller than M rec , the photon is combined with 
the corresponding fermion, i.e. the momenta of the photon and the fermion are added 
and associated with the momentum of the fermion, and the photon is discarded. 

3. Finally, all events are discarded in which one of the final-state fermions is within a 
cone of 10 degrees around the beams. No other cuts are applied. 

We consider the cases of a tight recombination cut M rec = 5 GeV and of a loose recom- 
bination cut M rec = 25 GeV. In the following observables, the momenta of the W bosons 
are always defined by the sum of the momenta of the two corresponding decay fermions 
after the eventual recombination with the photon. 

The results for the distributions have been obtained from 50 million events. In the 
following figures we always show on the left-hand side the absolute distributions in lowest 
order and including the corrections for the recombination cut M rec = 5 GeV, and on 
the right-hand side the corresponding relative corrections for the two recombination cuts 
M rec = 5 GeV and M rec = 25 GeV. 

In Figure [TJ we show the distribution of events in the angle between the W + and 
the incoming e + . Apart from the normalization effects, a distortion of the distribution 
occurs. This is mainly due to hard-photon emission from the initial state, which boosts 
the CM-system of the W bosons and causes a migration of events from regions with 
large cross section in the CM system to regions with small cross section in the laboratory 
system. This effect is also visible in the following distributions. The production-angle 
distribution hardly depends on the recombination scheme, as expected. An increase of 
the recombination cut leads to a small redistribution of events from the backward to the 
forward direction. 
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Figure 2: Decay-angle distribution for e + e — ► z/^/x+du and y/s = 200 GeV 



The distribution of events in the angle between the W + and the outgoing fi + is pre- 
sented in Figure |||. In this case we find a sizeable dependence on the recombination mass 
M rec for large decay angles, where the cross section is small. This originates from the fact 
that the recombination of a fermion with a photon parallel to this fermion decreases the 
angle between the fermion and the W boson from which the fermion results. 

The distribution of events in the energy of the outgoing /i + is depicted in Figure |^. 
In the on-shell approximation it would be restricted between 20.2 GeV < E^< 79.8 GeV. 
Outside this region, the corrections calculated in DPA are not reliable. The recombination 
of a photon with the muon increases the muon energy. Consequently, an increase of the 
recombination mass M rec shifts the distribution to larger muon energies as can be seen in 
the relative corrections. 

Finally, in Figures |] and |5] we show the distributions of events in the invariant masses 
of the final-state lepton pair, , and of the final-state quark pair, Md u . The depen- 
dence on the recombination cut is sizeable everywhere. The results for the invariant-mass 
distributions can be understood as follows. For small recombination cuts, in most of the 
events the W bosons are defined from the decay fermions only. If a photon is emitted 
from the decay fermions and not recombined, the invariant mass of the fermions is smaller 
than the one of the decaying W boson. This leads to an enhancement of the distribution 
for invariant masses below the W resonance. This effect becomes smaller with increas- 
ing recombination mass. The enhancement is proportional to the squared charges of the 
final-state fermion, i.e. it is largest for the leptonic invariant mass. On the other hand, if 
the recombination mass gets large, the probability increases that the recombined fermion 
momenta receive contributions from photons that are radiated during the W-production 
subprocess or from the decay fermions of the other W boson. This leads to positive cor- 
rections above the considered W resonance. The effect is larger for the hadronic invariant 
mass since in this case, two decay fermions (the two quarks) can be combined with the 



5 



dcr 



fb 
GeV 




40 60 
Ep [GeV] 



5 [% 

10 



-50 
100 



! 


1 1 1 

5 GeV - 


i 
i 

— i 

i 
i 
i 

i 
i 


j^^^ 25 GeV _ 

ft 
J' 

Tr 1^ 

i t 


i 

iL 

i 
i 
i 
i 


i Tl 
i i 

1 1i 


1 in' 1 






f \ 




i i 


i 1 



20 



40 60 
E^ [GeV] 



80 



Figure 3: Muon-energy distribution for e + e — > z^/i + du and ^Js = 200 GeV 
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Figure 4: Invariant-mass distribution of the lepton pair for e + e — > z/ A1/ u + du and \/s 
200 GeV 
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Figure 5: Invariant-mass distribution of the quark pair for e + e~ — > u fl fi + dn and y/s = 
200 GeV 



photon. The effect of the squared charges of the final-state fermions is marginal in this 
case because the contribution of initial-state fermions dominates. 

The distortion of the W invariant-mass distributions is of particular interest for the 
reconstruction of the W-boson mass Mw from the decay products. In order to illustrate 
the impact of the corrections on the determination of Mw, we fit the simple Breit-Wigner 
distribution 

/ da \ const. , . 

VdMV = (M 2 - M w>fit ) 2 + M Wifit r w>fit ( } 

to the W line shapes shown in Figures [| and [|, treating M Wjfit and r W fit (as well as the 
normalization constant) as free fit parameters. We determine the fitted W-boson masses 
from the predictions resulting from the lowest-order CC03 diagrams, M w °g^' CC03 , from the 
complete lowest-order diagrams, M w °g", and from the fully corrected predictions, M^ r J t . 
In addition we give the mass shift resulting from the corrections AM^j t = M^g t — M w °g°. 
The results of these fits, which are contained in Table |], show that the fitted W-boson mass 
changes at the order of some 10 MeV if the corrections are included. From the discussion 
of the line-shape distortion above it is clear that this mass shift is more positive if more 
photons are recombined. The results also illustrate that the fit results vary at the order 
of some 10 MeV for different fit ranges in M. 

We have also considered the case where each photon is recombined with the nearest 
fermion, i.e. M rec > -^i. In the distributions of the production angle, the decay angle, 
and the muon energy this increase of the recombination cut leads to the expected small 
changes of the corrections. For the invariant-mass distributions, the case of complete 
recombination differs from the case M rec = 25 GeV only by a slight decrease of the relative 
corrections. This is due to the fact that the recombination of a photon that forms an 
invariant mass of at least 25 GeV with charged fermions shifts the events by at least 
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M rec 


A# Born,CC03 












[GeV] 


[GeV] 


[GeV] 


[GeV] 


[MeV] 


W+ - 




5 


80.366 


80.365 


80.363 


-2 








80.372 


80.371 


80.365 


-5 






25 


80.366 


80.365 


80.374 


+9 








80.372 


80.371 


80.385 


+14 


W" 


— > du 


5 


80.365 


80.364 


80.378 


+ 14 








80.371 


80.370 


80.384 


+14 






25 


80.365 


80.364 


80.397 


+32 








80.371 


80.370 


80.415 


+46 



Table 4: Results for Mw.fit for a Breit-Wigner fit to the invariant-mass distributions 
shown in Figures [| and |^, using the fit ranges 78.3 GeV < M < 82.3 GeV (upper values) 
and 76.3 GeV < M < 84.3 GeV (lower values) 



3.8 GeV for the leptonically-decaying W boson and 7.4 GeV for the hadronically-decaying 
W boson. Such a shift leads to an increase of the corrections only outside the range 
shown in Figures § and [5[ Since there is no change in the shape of the distribution near 
the resonance, the results for the invariant-mass fit are essentially equivalent. 

In the plots on the right-hand side in Figures |] and |5| we have included the results for 
phase-space slicing and for the subtraction method. The agreement between these results 
demonstrates the correctness of our implementations and gives another estimate on the 
size of the integration errors. 

In summary, we have constructed the event generator RacoonWW for e + e _ — > 
WW — > 4/(+7) which includes the complete 0(a) electroweak corrections in double- 
pole approximation. With this generator we have calculated the total cross sections, 
including the CC03 cross sections, and various distributions of experimental relevance for 
typical LEP2 energies. The detailed numerical discussion of the corrections, in particular, 
illustrates the importance of the issue of photon recombination. More details of the 
RacoonWW approach and further results, including a comparison to results of other 
authors, will be presented elsewhere [ O . 
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